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We present a theoretical treatment of coherent light scattering from an interacting ID Bose 
gas at finite temperatures. We show how this can provide a nondestructive measurement of the 
atomic system states. The equihbrium states are determined by the temperature and interaction 
strength, and are characterized by the spatial density-density correlation function. We show how this 
correlation function is encoded in the angular distribution of the fluctuations of the scattered light 
intensity, thus providing a sensitive, quantitative probe of the density-density correlation function 
and therefore the quantum state of the gas. 

PACS numbers: pacs numbers 

The remarkable progress made in cold atom physics over the past two decades has been facilitated by the precision 
and flexibility of the measurement tools available to experimentalists. Optical probing by lasers has proved to be 
one of the most important tools, for example absorption and phase contrast imaging [l| have been extensively used 
to measure atomic density in bosonic and fermionic systems. Proposals for employing coherent light scattering from 
cold atoms are long standing 0]. More than a decade ago superradiant light scattering from a Bose condensate was 
demonstrated, and stimulated light scattering was used to implement Bragg scattering of matter waves More 
recently Mekhov et al Q have proposed light scattering as a probe of cold atoms in an optical lattice. 

In this letter we present a theoretical treatment of monochromatic light scattering appropriate for a ID Bose 
gas. We show that for far detuned incident light, a clear signature of various quantum states of the gas can be 
obtained from the fluctuations of the scattered light. The ID Bose gas is an important model in many-body physics, 
providing one of the simplest paradigms of a strongly interacting system, owing to its exact integrability via the 
Bethe ansatz Q. In addition, the distinction between fermions and bosons (that are well known in 3D systems), are 
suppressed in ID Q. For a ID Bose gas with repulsive interactions, the system state can be characterized by the two 
body spatial density-density correlation function Possible system states range from a nearly ideal gas in the 'high' 
temperature regime (where the spatial correlation function displays bunching), to the Tonks-Girardeau gas in the 
interaction dominated regime, (where the correlation function displays antibunching) . Experimental investigations of 
the ID Bose gas have been pursued in recent years A method for measuring momentum correlations of a cold 
atom system has been demonstrated by analysing the density fluctuations of a resonant absorption image [lo| . Two 
and three body correlations were also studied via resonant absorption images [oj . Impressive as these experiments are, 
this technique cannot image submicron density fluctuations. In this work we show that a nondestructive measurement 
of the spatial nonlocal density-density correlation function of a ID Bose gas can be made by measuring the angular 
distribution of the fluctuations in the scattered light intensity. 

We consider a system of two-state bosonic atoms in a long thin cylindrical trap, illuminated by a coherent laser 
beam of frequency wl and wavevector Icl. Each atom has resonant frequency ujo and transition dipole moment d. We 
use the rotating wave approximation for the radiation interaction. The laser mode is treated as a classical quantity 
Ej^g4(kL X--a;L)t^ and its interaction with the atoms is characterised by the detuning Al = uji^ — ujq and the Rabi 
frequency f^L = d • El/S. In the case where the detuning is large, Al ^ f^L,r, (P is the spontaneous emission rate) 
the upper internal atomic state can be adiabatically eliminated, leaving coupled equations for the scattered radiation 
modes a\^{t) and the atomic fleld operator for the lower internal state ^'(x). The equation for ai^{t) can be integrated, 
(ignoring small terms and employing the Markov approximation) leading to 

Mt) = ^^e— -*<5(Ak)e'^''*/^ / dh. e^^ ^pi^, t). (1) 

Here denotes the usual atom-radiation coupling constant, Ak = wl — Wk, p = ^^'5 is the 3D density of the system, 
q = Icl — k and S'^ is the region of space where the incoming light intersects the gas. We have excluded the additive 
quantum noise term from Eq. ([1]) as we will only consider normally ordered averages of a^. 

For the cylindrically symmetric trapping potential T^(x) = moj'^di'^ + z'^)/2, the ID regime can be achieved when 
fiujj_ ^ k^T,p, (the chemical potential). In addition the atomic recoil energy should be much less than huj^. In this 
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case the atomic field operator is well approximated by ^(x, i) = {^/nl±) 
4>{x,t) is the ID annihilation operator of the atomic ground state and Zj 
length. Performing the transverse integration in Eq. ([ij gives 



exp[— (y^ + z'^)/2l\ — iu±^t](}){x, t), where 
— ^h/(raijj±^ is the transverse oscillator 
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where A = 7rr2Lgke"^'''/^(5(Ak)e~(«y+''-)'-L/'^/AL, h = (j)^ (j) is the ID density and = [-L/2, L/2] is the region where 
the light intersects the gas. In this paper, we consider the case where the incident light is a plane wave polarized 
along the z axis, with Icl in the x-y plane and making an angle 0l to the y axis. The intensity of the scattered light 
is proportional to the total number of photons in a given mode, 



(«k«k) 
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N + 271^1 dr g'^'^\r) cos (q^r) {L - r) 
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where N = nL is the total number of atoms illuminated. The first term in the square brackets in Eq. ([3]) arises 
from the atomic shot noise, while the integral provides the effect of the atomic correlations, via the density-density 
correlation function of the gas, g'^-^^x — x') = (:n(x) n{x') We have assumed that, within the iUuminated region, 

the density is constant, {n{x)) = n, and g^^"^ is translationally invariant. 

The function g^^^ characterises the state of a uniform density ID Bose gas 0]. At equilibrium in the thermodynamic 
limit, ^'•^■'(r) is completely determined by two parameters [llj: the dimensionless interaction strength 7 ~ m2(jj±a/hn 
(where a is the 3D s-wavc scattering length, and m is the atomic mass) and the reduced temperature r = T/T^ where 
Td = h^v? / {2k^m) is the temperature of quantum degeneracy. Varying these two parameters gives rise to three 
physically distinct regimes, each of which can be further split into two subregimes. 

(i) Nearly ideal-gas regime, 7 ^ min{T^, -^/r}- Here the temperature always dominates over the interaction strength, 
and the local density-density correlation displays thermal bunching, 1 < g'^^\Q) < 2. This regime has two subregimes 
defined by r ^ 1 (quantum decoherent gas) or t ^ 1 (classical decoherent gas). 

(ii) Weakly interacting regime, <§; 7 ^ 1. This regime, for which both the interaction strength and the temper- 
ature are small, realises the quasicondensate where density fluctuations are strongly suppressed but long wavelength 
fluctuations of the phase still exist The density-density correlation is close to the uncorrelated value of g'^-'(r) w 1, 
and fluctuations occur due to either vacuum or thermal excitations. Two further subregimes can be defined, r ^ 7 
(dominant vacuum fluctuations) or r ^ 7 (dominant thermal fluctuations). 

(iii) Strongly interacting regime (Tonks gas), 7 3> max{l,\/T}. Here the interaction strength always dominates 
over temperature induced effects, and the local density-density correlation displays interaction-induced antibunching 
< 5*'^''(0) < 1- This regime is usefully divided into two further subregimes defined by r 1 (low temperature Tonks 
gas) and r ^ 1 (high temperature Tonks gas). 

In each of the above regimes, perturbation theory can be applied and six different analytic expressions for the density- 
density correlation function have been calculated 



ll| . The density-density correlation function decays from its local 



value g''^'(0) to the uncorrelated value g'-^' = 1 over a length 1^ ( the correlation length). For the quantum decoherent 
gas we have Ic — 2/(™) (the phase coherence length), while for the classical decoherent gas, Ic = \J Att / [rn^) (the 
thermal de Broglie wavelength). For the weakly interacting quasicondensate regime, (with either vacuum or thermal 



fiuctuations dominant) we have Ic 



Ic ^ l/n and finally for the high temperature Tonks gas we have I. 
The physically relevant cases for light scattering are where L ^ I 



(the healing length). For the low temperature Tonks gas we have 



the coherent scattering it is convenient to define ri{r) = g^'^\r) — 1, Eq. ([3]) then becomes 



. . . 

> A. In order to distinguish the incoherent from 



(4«k) = (4)(«k) + i^i'^^ 
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dr Tj (r) cos (q^r) {L 
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The first term on the right of Eq. Q , (a]^) (oj,) , is the coherent scattering and is predominantly within angular width 
^ \/ L oi the (jx = direction. For X L, the total intensity of coherently scattered light at a distance r from the 
gas is 



Ic = N' 



3sin''6' 
' 87rr2 



-smc 



(5) 
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where 6 = k^j. — /cl sin 6*008 6*5, and 9^ = arctan(|^) is the in-plane scattering angle, and 9 ~ arctan(^^^-| — -) is the 
out-of-plane scattering angle (see Fig. 2]). 

The signature for the atomic g^^\r) function is in the fluctuations of the scattered intensity arising from the integral 
term of Eq. (|4]). These fluctuations are typically expressed in terms of quadratures of the radiation field operator 

h.c.) . The variances of these 



(relative to a particular phase reference 
quadrature operators. 



), which are defined by Xg 



\A\'^N 



N 
I? 



t; + T2 '?(^) cos{q^r){L-r) + cos(2/3) 



sin [q^{L - r)] 
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can be measured by well known techniques. For L ^ Ic > X, the dependence on /3 in Eq. 1^ is very tightly concentrated 
around the transmitted and reflected directions; i.e. \qx\ ^ 1/L. A simple and prudent choice of phase reference is 
(3 = J. In this case the quadrature noise simplifies to w jApiV 1/2 + n J^dr ri{r) coslq^r) , which consists of the 

atomic shot noise term, and a Fourier transform of r]{r). 

In order to obtain an analytic understanding of our results, we first take the limiting cases where 7 and t take on 
values which lie deep within one of the regimes defined above. The density-density correlation then simplifies to 



?7(r) 



exp [— r?ir] (a) 

exp (b) 

— exp [— r7i^r^/2] (c) 

— sin^(7r??r)/(7rnr)^ (d) 



(a) is deep within the decoherent quantum regime, (b) is deep within the decoherent classical regime, (c) is deep 
within the high temperature Tonks regime, (d) is deep within the low temperature Tonks regime. The integral for 
R\ can then be calculated analytically; 



l + 2n^T/{ql + n'^T'^) 
\A\^N J 1 + 2v A^7^ e~'?x/(2»^r) 

1 - 27^^7^6-9^/(2"'-) 
1 -tri [gx/(27rn)] 



(a) 
(b) 
(c) 
(d) 
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where tri(7-) ~ max(l — |?'|,0). We note that for = 0, the case of a perfectly coherent quasicondensate, density 
fluctuations are completely suppressed, and R^ reduces to the shot noise. For intermediate values of 7 and r, we 

solve for R^ numerically using expressions for g^^) taken from [llj. As seen in Fig. 21 the intensity fluctuations 
arising from ri{r) arc typically scattered over a wider angle than the coherent scattering. For example, the thermal 
case in Eq. ©(b), the angular spread is ~ arctan (Af(A/L)(2.3-y/r/27r)). As a crude approximation, we can write 
r]{r) = rje^^ where —1 < ?? < 1, and we find the incoherent scattering intensity, at a distance r from the gas, to 

be lir = 



N 



2At 



Thuj\ 



3 sin^6 
' 87rr2 



1 + ?i^c?7v27r expl 
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Intensity isosurfaces of both the coherent and incoherent scattering are shown in Fig. 2] 

The results for i?^ plotted against the scattered angle are shown in the right-hand columns of Figs.[TH3l for a range 
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of T and 7. In the left hand-columns we plot the corresponding g^'^\r). In Fig. [T] the sequence of plots are at a fixed 
high temperature, and show the transition from a decoherent classical gas (top row) to a high temperature Tonks 
gas, as the interaction strength 7 increases from to 10. For weak interaction g'^) displays a peak at r w 0; thermal 
bunching. As 7 increases, the peak in y'-^) changes to a dip; antibunching. The behaviour of the corresponding 
scattered intensity fluctuations, which consist of a noise floor {\A\'^N/2)^ plus the Fourier transform of ry = g^^) _ 

can be readily understood. A positive (negative) function ?/ produces a peak (trough) in i?^ , with a width scaling 
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inversely to Ic- This is clearly illustrated in the sequence of plots for 0l = in Fig. [TJ Atom bunching produces 
enhanced intensity fluctuations over a range of angles around 0^ = 0, while antibunching results in noise suppression 
in the scattered light over a similar angular range. We can interpret the more complex structure in Fig 1 (c2) as 
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FIG. 1: (color online). Density-density correlation function (LH column) and corresponding scattered intensity quadrature 
noise (RH column) . Sequence shows transition from a decoherent classical gas to a high temperature Tonks gas. Incident angle 
= (solid blue lines); 7r/4 (red dashed lines); 7r/2 (green dot-dashed lines). Other parameters (kL)x = 20n, and L 3> Ic- 
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FIG. 2: (color online). Density-density correlation function and scattered intensity quadrature noise for the transition from a 
decoherent classical gas to a zero temperature quasicondensate. Incident angle = (solid blue lines); = 7r/4 (red dashed 
lines); = -k 12 (green dot-dashed lines). (kL)x = 5n, and L 3> Ic- 



being due to a remnant bunching on a long spatial scale superimposed on a more dominant (and spatially narrower) 
antibunching. For an incident angle of 7r/4, these features are replicated around a (broader) transmitted peak at 
= — 7r/4 and a reflected peak at — 37r/4. Fig. [21 where all plots are at weak interaction strength 7 = 0.01, shows 
the transition from a decoherent classical gas to a zero temperature quasicondensate as temperature is decreased to 
zero. At the higher temperatures, the scattered intensity fluctuations are enhanced in a broad angular region around 
the transmitted (and reflected) beams. At the lowest temperature, where g^"^^ displays weak antibunching with 1^ 
given by the condensate healing length (which is long), the light fluctuations arc suppressed over a narrow angular 
range. Fig. [^^bl) exhibits competition between bunching and antibunching on different length scales, resulting in a 
subtle minima in g^^' at nr w r/(47). This leads to a sharp dip in the noise features Fig. [Ub2) showing that light 
scattering is an ideal probe of these subtle, and controversial [ill, density-density correlations. Finally, in Fig. |3l a 
sequence of plots at constant large interaction strength 7 = 8, shows the transition from a high temperature to low 
temperature Tonks gas. In our analysis we have assumed the atoms are in an equilibrium state unperturbed by the 
interaction with the light. This assumption is valid provided the heating per atom from the recoil energy is much less 
than the initial energy per atom, Eq. The recoil component perpendicular to the trap axis is absorbed by the trap as 
a whole (because it is very tight), and so for 6'l = 0, only the incoherently scattered light contributes to the heating. 
We can estimate the heating rate per atom to be T?i^(fca;)^/2m where (fca;)^ is the mean square recoil momentum and 



T = r(^^)^ is the rate of photon scattering [131 . The perturbation to the initial state of the gas will be negligible 
provided Ti <C £'o- Expressions for Eq can be found from 
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FIG. 3: (color online). Density-density correlation function and scattered intensity quadrature noise for the transition from a 
high temperature Tonks gas to a low temperature Tonks gas. Incident angle = (solid blue lines); = 7r/4 (red dashed lines); 
= 7r/2 (green dot-dashed lines). (kL)x ~ 20n, L ^ Ic- 




FIG. 4: Intensity isosurfaces for coherently and incoherently scattered light. In each of the above figures, the blue cylinder 
shows the orientation of the ID Bose gas, the red cylinder shows the direction of the incoming laser beam, and the semi- 
opaque, purple surface shows an intensity isosurface as a function of scattering angle. In (a) we show the intensity of the 

coherently scattered light [given by Eq. (5) in the main text] with A'^^ (^^^) Thwi^ ^^^.j^ = 1, /clx = 0.7A;l and ki^L ~ 28. In 

(b)-(e) we show the incoherent intensity [given by Eq. (8) in the main text], for the approximation ri(r) = -qe ^ with 

^ (^r) rficJL gj^^y ~ 1, fcLx ~ 0.7kh, nlc = 1 and khlc ~ 7. In (b) we have rj — 0.5; in (c), 77 — 0.3; in (d), rj — —0.3; and in 

(e), J7 = -l. 



In summary, we have shown that the quantum state of a ID Bose gas can be nondestructively measured using off 
resonant light scattering. The angular distribution of the quadrature noise in the scattered light provides a direct 
signature of the nonlocal two-body correlations which characterize different system states of the gas. The method 
is useful across the full range of temperatures and atomic interaction strengths, and is capable of revealing subtle 
features in the density fluctuations that occur across multiple length scales. In addition, the method could also be 
used to monitor the dynamics of two-body correlations in the system, by simply solving the inverse Fourier transform 
of the intensity fluctuations. 
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